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Below, we detail the calculation of the correlation functions for the mem- 
O brane composition and deformation appearing in the main text (Eqs.(5,6)). 

We first define the integral: 
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which has the asymptotic behaviour 
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We are interested by the correlation function of a variable satisfying a 
dynamical equation of the form (Eq.3, main text) 
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yielding the correlation function 
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from which the results given in Eqs(5,6) can easily be derived. 

The membrane response to a point force being radially symmetric, the 
real-space correlation functions (such as those shown in Fig. 2 of the text) 
can be calculated with 
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where J is the Bessel function of the first king, satisfying J Q 27r d6e iq ' rcosd = 
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Figure 1: Reproduction of the phase diagramm appearing in the text (Fig.l) 
with indication of the states represented in the different supplementary 
movies (green stars, number corresponds to the movie numbers). 
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Microphase separation of membrane components is thought to play an important role in many 
physiological processes, from cell signaling to endocytosis and cellular trafficking. Here, we study 
how variations in the membrane composition can be driven by fluctuating forces. We show that the 
membrane steady state is not only controlled by the strength of the forces and how they couple to 
the membrane, but also by their dynamics: In a simple class of models this is captured by a single 
a correlation time. We conclude that the coupling of membrane composition to normal mechanical 
forces, such as might be exerted by polymerizing cytoskeleton filaments, could play an important 
role in controlling the steady state of a cell membrane that exhibits transient microphase separation 
on lengthscales in the 10-100 nm regime. 



Biological membranes are composed of a large number 
of different lipid and protein species, which are known 
to phase separate into macrodomains when isolated and 
cooled below a critical temperature [1-3]. The existence 
of nano scaled membrane domains at physiological tem- 
perature (including so-called lipid rafts) have been pos- 
tulated to play an important role in the clustering of sig- 
naling proteins [4] . The lateral organization of biomem- 
branes has in particular been linked to the structure of 
the cell cytoskeleton [5-7], and to membrane trafficking. 
It is now becoming increasingly clear that the control- 
ling principles of biomembrane organisation include the 
stochastic dynamics of its environment, which influences 
the membrane through biochemical signaling, mechan- 
ical perturbation, and direct material fluxes. Recently 
we showed how recycling of membrane components can 
sensitively tune the state of the membrane [8]. Here, 
we analyze the response of the membrane to a fluctuat- 
ing environment, and show that the level of membrane 
organization can be tubed by correlations in the dynam- 
ical environment. We first calculate the response to a 
locally correlated stochastic coupling of arbitrary origin, 
and show how the correlation time of the perturbation 
controls the membrane's response. We also analyze in 
detail the case in which membrane composition is cou- 
pled to the local membrane curvature, itself driven by 
random normal forces generated by the cytoskeleton. 

The most physiologically relevant thermodynamic 
transition in biomembranes is thought to be the liquid 
disordered (LD) to lipid ordered (LO) transition, which 
typically separate saturated and unsaturated lipids, and 
is promoted by cholesterol [4]. Since the LO phase ex- 
hibits reduced lipid mobility compared to the LD phase, 
in-vivo studies of the organization of biological mem- 
branes can emploit an approach based on dynamical 
tracking of tagged membrane component [5]. A par- 
ticularly interesting example [9] demonstrates the exis- 
tence of short-lived (~ 250msec.) low mobility nano- 
domains (<~ 75nm), appearing and disappearing near 



seemingly well-defined positions on the membrane. Such 
behaviour could result from biochemical or physically lo- 
calized stochastic forces that act to drive changes in mem- 
brane composition. 

Model for membrane phase separation. Our start- 
ing point is to write down a Hamiltonian in terms of the 
membrane composition field p(r) and normal membrane 
displacement u(r), both assumed to represent small per- 
turbations around a perfectly mixed (p = 0) and flat 
(u = 0) membrane [10]: 

+n(V 2 u) 2 + cr(Vu) 2 - 2kC' oP V 2 u - 2(p - 2fu) (1) 

which can either be motivated on symmetry grounds 
or, by tracing the physical interpretation of the various 
terms. These include: a local mixing interaction and 
concentration gradient term for the membrane composi- 
tion (parameters b and p), the membrane tension a and 
bending rigidity k, and a coupling term C' Q capturing 
how the local membrane spontaneous curvature depends 
on the composition to lowest (linear) order in p. Thus 
p is the density of membrane component (s) that couple 
to curvature. The fields / and ( are used to capture the 
coupling of external (cellular) forces conjugate to u and p 
respectively. In the following, we will mostly work in the 
spatial Fourier space, with x q = J c? 2 re lq r x(r), giving 
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The kinetic evolution of the two fields, subjected to 
local fluctuating forces /(f) and ((t) coupling to the dis- 
placement and composition, respectively, is given by the 
equations [11]: 

p q + Aq 2 (h qPq - /3 q u q ) = Aq 2 ( q 
r]qu q + (k q u q - (3 qPq ) = f q (3) 
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where A is related to the diffusion coefficient of membrane 
component by D = Ah q (D Q = D(q — »• 0)), and where r] is 
the solvent viscosity dampening the membrane displace- 
ment [12]. Our goal is to study the effect of correlated 
fluctuating forces (of zero mean) on the membrane state, 
to which end we adopt 



disregard the coupling between membrane composition 
and deformation (f3 — 0) and study the fluctuations of 
membrane composition subjected to a correlated fluc- 
tuating force £(t). The solution of Eq.(3) subjected to 
Eq.(4) can then be calculated (sec Supplementary Infor- 
mations - S.I.): 
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In the case of more than one localised force we can write 
/(r,t) — J2ifi( r ~ , t) , with fi the force due to the 
cytoskeletal anchor; similarly for £ where its canon- 
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ical force would be due to some localised biochemical 
coupling(s). Thus, with /,(r - r;,£) = fi(t)S(r - 
(here we consider applied forces / with zero range spa- 

tial correlations) f q _= e*^ /j (i) with (/<(«)) = and satisfying 'b q \ V 1 are" essentially unaffected, while 



This expression gives the appropriate static response 
{Pq = Cq/hq) in the limit Dq 2 r >• 1, and satisfies equipar- 
tition of energy: (|p g | 2 ) = k B T/h q in the limit of tem- 
porally uncorrelated forces (r p — ► 0) if the fluctuating 
force satisfies the fluctuation-dissipation theorem: 
(IC 9 | 2 ) = Co = K T /( A( l 2 Tp)- Eq.(5) shows that modes 
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stant. In what follows we restrict our attention to single 
force centres situated at r — noting that the membrane 
response is linear and so membrane deformation caused 
by additional force centre(s) would simply add indepen- 
dently. For a numerical solution of the microphase sep- 
aration induced by multiple force centres see supporting 
movies 1-4 while the level of interactions between these 
force centres is shown in Fig.l. 

Fluctuations of membrane composition. We first 
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FIG. 1. The dynamics of a membrane containing multiple 
force centres (of density p) depends on the ratio of the mem- 
brane length A (see text) to the distance between force cen- 
tres: Xy/p, and the ratio of the composition diffusion length 

As we 



V Dt to the distance between force centres ^/Drp. 
discuss in the text independent force centres each contribute 
linearly to the membrane displacement and enrichment. Each 
cartoon represents two neighbouring force centres, with a typ- 
ical separation p~ 1,/2 , and shows the diffusion length in red, 
the membrane length in blue (not to scale) and the region of 
substantial compositional variation in dark grey shading. In 
each case the controlling lengthscale is shown as a continu- 
ous coloured line and the other lengthscale is shown dotted. 
The shaded regimes shows where domains are substantially 
non-interacting. 



the fluctuating nature of the correlated force. 
Membrane-cytoskeleton coupling. The cytoskeleton 
underlying the plasma membrane has long been thought 
to influence the mobility and clustering of membrane pro- 
teins. So far, attention has mostly been directed to- 
ward the role of cytoskeleton filament lying parallel to 
the membrane, thereby creating hindrance to the lateral 
motion of membrane proteins (the so-called picket-fence 
model) [13]. The membrane composition is however also 
known to be sensitive to the local membrane conforma- 
tion (mostly its curvature [14]), which should fluctuate 
strongly under the action of mechanical forces produced 
by polymerizing and depolymerizing cytoskeleton fila- 
ments. This situation can be mathematically addressed 
through the parameter (3 coupling the local membrane 
composition to its local curvature. One can in principle 
solve Eq.(3) for the two coupled fields of membrane com- 
position and deformation (see for instance [11]). For the 
sake of clarity here we only investigate the limit of rather 
weak coupling (/3 2 < k q h q ), so that the membrane defor- 
mation is substantially determined by the external force, 
with little feedback from the membrane composition. It 
is known that, in the opposite limit of strong coupling, a 
curvature instability can arise where large membrane un- 
dulations coupled to strong composition variation grow 
spontaneously [10]. 

The Fourier transform of the curvature of a membrane 
deformed by a static point force is q 2 u q = f/(a(l + 
(Xq) 2 )), with A = \] njo the characteristic "healing 
length" of membrane deformation (~ lOOnrn for k = 
20fc B T and a = 10~ 5 J/m 2 ). In real space, the curvature 
is proportional to the modified Bessel function of the sec- 
ond kind K (r/X) [15], decaying exponentially at large 
distance, but diverging logarithmically at short distance. 
Solutions must be regularized by introducing a short- 
length cut-off a, effectively corresponding to the (molec- 
ular) lateral size over which the point force acts (we as- 
sume a <C A) . This linear description is valid in the small 
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FIG. 2. a Comparizon of the membrane dynamical response (as defined by Eq.(8)) to a fluctuating environment of typical 
correlation time r that drives the membrane composition changes via (i) localized exponentially decreasing forces (s cxp ), 
(ii) localized forces with guassian spatial distribution (s gauss ) and (iii) a local coupling between membrane composition and 
membrane curvature, itself excited by localized normal mechanical forces (s cur v)- The inset shows the different behaviors of 
the response functions for weak correlations: r — > 0. b. Local correlation function (p(r) 2 } driven by a fluctuating random 
force with a correlation time r (shown for r = 0.01, 0.1, 1, 5, 10). Left shows the case of a local mechanical driving force 
with coupling between membrane composition and curvature, Right: shows the case of a force decaying exponentially from 
a local source. In dashed is shown the expected result for an infinitely correlated (permanent) force, c. Evolution of the 
amplitude at the origin (p(0) 2 ) (black) and the distance ri/2 from the origin where the amplitude has dropped by half (grey) 
for an exponential perturbation, as a function of the perturbation correlation time r. The global measure C (dashed grey) is 
numerically shown to follow the expected scaling behaviour C = (p(0) 2 }7- 2 , 2 (dots) 



deformation limit Vu <C 1, requiring that the strength of 
the point force is smaller than y/na ~ 10 — 20pN, beyond 
which the membrane deforms into a tubular spicule. 

The mean squared curvature solution of Eq.(3) under 
the fluctuating force / satisfying Eq.(4) is (see S.I.): 
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where a q = k q /{riq) is the relaxation rate of the mem- 
brane deformation [16]. As a result of its coupling to the 
local curvature, the membrane composition sees an addi- 
tional fluctuating force with correlation ((' q (t)C' q (t+5t)) = 
Pq (u q (t)u q (t + St)} coming from the fluctuations of the 
membrane shape. 

The composition correlation function now involves 
three relaxation rates: of the force correlation 1/r, of the 
membrane deformation a q , and of the membrane compo- 
sition Dq 2 . In the regime of most interest to us the mem- 
brane shape adjusts much faster than the membrane com- 
position a q 3> Dq 2 , which is the case provided the mem- 
brane length A is much smaller than n/(Dri)(^ 10/im). 
This is always the case in practice (A ~ 10~ 2 — 1/im de- 
pending on membrane tension), the membrane deforma- 
tion can therefore be considered to relax instantaneously 
to its mechanical equilibrium on the timescale of compo- 
sition changes. The correlation functions then reads: 
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akin to Eq.(5) but with an effective driving force that 
includes the membrane elasticity and the strength C 



of the composition coupling to curvature but which is 
insensitive to any details regarding the dynamics of its 
deformation. 

Quantification of membrane organization. The fulll 
spatio-temporal correlation function for the membrane 
composition (see S.I.) can, in principle, yield any statis- 
tical information of the membrane composition, e.g. by 
inverse transformation to real space (Fig. 2b). This step 
typically involves numerical calculation so it is desire- 
able to also obtain an analytic statistical measure for the 
membrane state. Towards this end we define the inte- 
grated autocorrelation function: 
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here normalized by the membrane's response to a per- 
manent perturbation (C(r p — > oo) = 1) to reveal the role 
of the fluctuating nature of the force. In principle, the 
susceptibility of the membrane composition h q = b + jiq 2 
also contains the static composition correlation length 
\J fi/b. We consider the membrane to be far from the 
critical point for membrane composition, in which case 
this length scale is of molecular dimension [11] and we 
may use the approximation /i = henceforth. 

Eq.(8) can readily be used to quantify the membrane 
composition's response to fluctuations of the local curva- 
ture under stochastic mechanical forces (using Eq.(7)): 



C — 



t(t — log T — 1) 



(9) 



with t = Dt/X 2 . This results can be compared to the 
response of other types of fluctuating environments, such 
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as an exponentially decreasing local source of signaling 
molecules (with ( q = 1/(1 + (Aq) 2 ) 3 / 2 ): 

_ f(f 2 -2flogf-l) 

or a similar force of Gaussian form: ( q = e^^ q ^ / 2 . The 
corresponding response functions are shown in Fig. 2a. 
All three show a slow convergence toward the static re- 
sponse (C — > 1 for f — > oo), but exhibit quite distinct 
behaviour for weakly correlated excitations (small f): 

^curv ~ ^ 7" log — C C xp ^ ^gauss ~ ^ T ( 1 1 ) 

r— >0 T T— >0 

This difference comes from the fact that the membrane's 
conformational response to a localized mechanical force 
contains a high density of short wavelength modes to 
which the membrane composition is quick to adjust by 
diffusion, thereby leading to a logarithmically diverging 
sensitivity to correlations of the stochastic environment 
for C curv . Note that Eqs.(9,10,ll) are valid provided r > 
a 2 /D, which is a weak constraint since the short-distance 
cutoff a is of molecular size. 

Two important measures of the fluctuating membrane 
"domains" are the intensity of phase segregation within 
them, and their average size. Both properties depends 
on the correlation time r, and are somehow mixed in the 
global measure C of Eq.(8). As show in Fig. 2b, the real- 
space correlation function (p(r,t) 2 } can be obtained from 
the (numerical) inverse Fourier transform of (p q (t)p q >(t)} 
(see S.I., Eqs.(5-6)). These plots show that while the 
characteristic shape of the domain resembles the perma- 
nent shape even for short correlation time r (of order 
~ X 2 /D), the amplitude of the reorganization requires 
longer times to reach its static level. This behaviour 
is further illustrated in Fig. 2c, where both the maximal 
level of correlation {(\p{a)\ 2 )) and the typical distance 
r\j2 from the source at which the amplitude falls by half 
((|y°( r i/2)| 2 ) = (|p( a )| 2 )/2), are plotted as a function of 
the correlation time. Remarkably, the global measure of 
the perturbation C closely follows the expected scaling re- 
sult C <~ (p(a) 2 )ry 2 for all times, thereby validating C as 
a good quantity with which to measure the membrane's 
response to a fluctuating environment [17]. 
Conclusion. The processes that control the lateral orga- 
nization of biological membranes are still far from being 
well understood. Much effort has been put into under- 
standing the equilibrium phase behaviour of these mem- 
branes as part of a philosophy that sees them as quasi- 
equilibrium structures. We believe that active fluctua- 
tions in the membrane's environment play a crucial role 
in this microphase separation and that this facilitates 
the sensitive transmembrane signalling of environmental 



variations. Our model shows the crucial importance of 
the temporal correlation of a noisy environment in trig- 
gering the membrane response. Fig.l 

Fig. 2 shows the relative strength of the response with 
the correlation time of the perturbation (compared an 
infinite correlation time). While the extent of the mem- 
brane reorganization reaches the extent of the driving 
force even in a highly fluctuating environment (for cor- 
relation times larger than X 2 /D, of order 100 msec for 
A = lOOnm and D = 0.1/xm 2 /s), the amplitude of the 
perturbation must be driven much more slowly if it is 
to reach its full potential. A primary result of this work 
is that compositional changes due to variations in mem- 
brane curvature driven by normal forces exerted by the 
cytoskeleton could represent the dominant cellular strat- 
egy for membrane organization. 
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